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We review cosmology in non-minimal Yang-Mills/Maxwell theory, in which the Yang- 
Mills/electromagnetic field couples to a function of the scalar curvature. We show that power- 
law inflation can be realized due to the non-minimal gravitational coupling of the Yang-Mills field 
which may be caused by quantum corrections. Moreover, we study non-minimal vector model in 
the framework of modified gravity and demonstrate that both inflation and the late-time acceler- 
ated expansion of the universe can be realized. We also discuss the cosmological reconstruction of 
the Yang-Mills theory. Furthermore, we investigate late-time cosmology in non-minimal Maxwell- 
Einstein theory. We explore the forms of the non-minimal gravitational coupling which generate 
00 . the finite-time future singularities and the general conditions for this coupling in order that the 

' finite-time future singularities cannot appear. 

(N 

> ■ I. INTRODUCTION 

o 

I^H I It is observationally confirmed that the current expansion of the universe is accelerating, as well as that there 
I existed the inflationary stage in the early universe [J, 01- The former phenomenon is called "dark energy" problem 
(for recent reviews, see [1, Ij])- The studies of the dark energy problem can be categorized into the following two 
approaches. One is to introduce some (unknown) matter which is responsible for dark energy in the framework of 

, general relativity. Another is to modify the gravitational theory, e.g., to study the action described by an arbitrary 

^ ' function F{R) of the scalar curvature R, which is called "F(i?) gravity" (for a review, see 0, Such a modified 
^ ' gravity must pass cosmological bounds and solar system tests because it is considered as an alternative gravitational 
(~| [ theory. 

A very realistic modified gravitational theory which evades solar-system tests has recently been proposed by Hu 
' and Sawicki (for related studies, sec Q). In this theory, an effective epoch described by the cold dark matter 
^ , model with cosmological constant (ACDM), which accounts for high-precision observational data, is realized as in 
' general relativity with cosmological constant. Although this theory can successfully explain the late-time acceleration 
, of the universe, the possibility of the realization of inflation has not been discussed. In Refs. [1, H, [I^l, following the 
T-H ■ previous inflation-acceleration unification proposal modified gravities in which both inflation and the late-time 
acceleration of the universe can occur have been investigated. The classification of viable F{R) gravities has also been 
suggested in Ref. [§]. 

As another gravitational source of infiation and the late-time acceleration of the universe, there is a coupling between 
OO . the scalar curvature and matter Lagrangian [13, (see also [13 )■ Such a coupling may be applied for the realization 
' of the dynamical cancellation of cosmological constant [1^1 . In Refs. [igI. [TtI. [isj . the criteria for the viability of such 
theories have been examined. As a simple case, a coupling between a function of the scalar curvature and the kinetic 
term of a massless scalar field in a viable modified gravity has been considered in Ref. [l^ . 

It is known that the coupling between the scalar curvature and the Lagrangian of the electromagnetic field arises in 
' curved spacetime due to one-loop vacuum-polarization effects in Quantum Electrodynamics (QED) [2^. In the present 
article, following the considerations in Ref. [2l|, we review cosmology in non-minimal non-Abelian gauge theory (Yang- 
Mills (YM) theory), in which the non-Abelian gauge field (the YM field) couples to a function of the scalar curvature. 
In particular, it is shown that power-law infiation can be realized due to the non-minimal gravitational coupling of the 
Yang-Mills field. The consequences presented correspond to the generalization of the results for non-minimal Maxwell 
theory with the coupling of the electromagnetic field to a function of the scalar curvature [l^. Furthermore, we 
consider a non-minimal vector model in the framework of modified gravity. It is demonstrated that both inflation and 
the late-time accelerated expansion of the universe can be realized. Wc also study the reconstruction of the YM theory. 
In the past studies, inflation driven by a vector filed [1^ and its instability [2^ and gravitational-electromagnetic 
infiation from a 5-dimensional vacuum state [2^ have been discussed. As a candidate for dark energy, the effective 
YM condensate [1^, the Born-Infeld quantum condensate [l^ and a vector field [l^l have been proposed. Models of 
vector curvaton have also been constructed [1^ . 

In addition, following the investigations in Ref. (soj . we review late-time cosmology in the non-minimal Maxwell- 
Einstein theory. Wc investigate the forms of the non-minimal gravitational coupling of the electromagnetic flcld 
generating the finite-time future singularities and the general conditions for the non-minimal gravitational coupling 
of the electromagnetic field in order that the finite-time future singularities cannot appear. In Ref. [3l|, F{R) gravity 
coupled to non-linear electrodynamics has been examined. 
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This article is organized as follows. In Sec. II we investigate a non-minimal gravitational coupling of the SU{N) 
YM field in general relativity. In Sec. Ill we consider non-minimal vector model in the framework of modified gravity. 
In Sec. IV we discuss the reconstruction of the YM theory. In Sec. V wc study non-minimal Maxwell-Einstein theory 
with a general gravitational coupling. We explore the cosmological efi'ects of the non-minimal gravitational coupling of 
the electromagnetic field to a function of the scalar curvature on the finite-time future singularities. Finally, summary 
is given in Sec. VI. We use units in which kB = c = h = 1 and denote the gravitational constant SttG by k^, so that 

w 

adopt Heaviside-Lorentz units. 



k2 = 87r/Afpi^ where Mpi = Q-^/^ = 1.2 x IQi^GeV is the Planck mass. Moreover, in terms of electromagnetism we 



II. INFLATION IN GENERAL RELATIVITY 



In this section, we study non-minimal YM theory in general relativity. The model action is given as follows 21|: 

S'gr = [ d'^Xy/^{CEii+ Cym) , (2.1) 

(2.2) 

Cym = --/(i?)f,?„F°^" 1 + bg'ln ^ ^ ^ , (2.3) 

m = l + f{R), b^^^^N, F^.^d.Al^d^Al + r'-A^^Al, (2.4) 

where g is the determinant of the metric tensor p^i,, R is the scalar curvature arising from the spacetime metric tensor 
g^v, and >Ceh is the Einstein-Hilbert action. Moreover, £ym with I{R) = 1 is the effective Lagrangian of the SU{N) 
YM theory up to one-loop order [H, [s^, f{R) is an arbitrary function of R, b is the asymptotic freedom constant, 
F^j, is the field strength tensor, A° is the SU{N) YM field with the internal symmetry index a (Roman indices, a, 
b, c, run over 1,2, .. . — 1, and in F^^^F"'^'^ the summation in terms of the index a is also made), and /"'"^ is a 
set of numbers called structure constants and completely antisymmetric [sj]- Furthermore, /i is the mass scale of the 
renormalization point, and a field-strength-dependent running coupling constant is given by (ssj 



1 + 6.g2 In 



X = — F^'^F"'"' . 



(2.5) 



where g is the value of the running coupling constant when X = fi^. 

Taking the variations of the action (|2.1[) with respect to the metric g^^, and the SU{N) YM field A° , we obtain the 
gravitational field equation and the equation of motion of as 
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and 
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{y^I{R)eF''^"') - I{R)ef'''"'AlF''t"' = , 



(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2.10) 



respectively, where T^^^' is the contribution to the energy-momentum tensor from the SU {N) YM field, the prime 
denotes differentiation with respect to R, is the covariant derivative operator associated with g^^, □ = .g^'^V^Viy 
is the covariant d'Alembertian for a scalar field, and Rfj_i, is the Ricci curvature tensor. Moreover, e is a field-strength- 
dependent effective dielectric constant [33| . and e « 2.72 is the Napierian number. In deriving the second equalities 
in Eqs. and (EH), we have used X = - (1/2) F/^^F^^''. 

We assume the flat Friedmann- Robertson- Walker (FRW) spacetime with the metric ds^ = ~dt^ + a^{t)dx^ = 
a^(r]){—drf' + dx^), where a is the scale factor and -q is the conformal time. 
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In the search of exact solutions for non-minimal YM (electromagnetic)-gravity theory, the problem of off-diagonal 
components of YM (electromagnetic) stress tensor being non-zero while the right-hand side of Einstein equations is 
zero (for the argument about the problem of off-diagonal components of electromagnetic energy-momentum tensor 
in non- minimal Maxwell-gravity theory, see (ssj). As a simple case, we can consider the following case in which the 

off-diagonal components of T^^^' in Eq. (|2.7p vanishes: (i) Only (YM) magnetic fields exist and hence (YM) electric 
fields are negligible, (ii) JB" = {Bf, , ), where Bf = i?f = 0, ^ 0, namely, we consider the case in which only 
one component of JB° is non-zero and other two components are zero. In such a case, it follows from divB" = that 
the off-diagonal components of the last term on the right-hand side of T^^^', i.e., V^Vi^ {f'{R)!F) are zero. Thus, all 

of the off-diagonal components of t'i^^^ are zero. Throughout this article, we consider the above case. 

In Eq. (|2.10p , we can neglect the higher order than or equal to the quadratic terms in A° because the amplitude of 
is small, and investigate the linearized equation of Eq. (|2.10p in terms of A° in the Coulomb gauge, d^A^{t, x) = 0, and 
the case of Ag(i, x) = 0. Replacing the independent variable t by rj, we find that the Fourier mode Af{k, rj) satisfies 
the equation {d^Af{k,r])) / {drf-) + {1/ I{r])) [dl{r])/ {dr])][dAf{k,r])/ [dr])] + k^Afik,r]) = 0. By using the WKB 
approximation on subhorizon scales and the long- wavelength approximation on super horizon scales, and matching 
these solutions at the horizon crossing [36| . we find that an approximate solution is given by 



\AUKv)\' = \c{k)\' = 



2kl{r,k 



dliVk) 



2fc/(77fc) dr, 



+ i]k 
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(2.11) 



where rjk and rjf arc the conformal time at the horizon-crossing and at the end of inflation, respectively. It follows 
from Bi^'^°^°'^ {t,x) = a~^eijkdjAk{t,x), where e^fe is the totally antisymmetric tensor (ei23 = 1), that the amplitude 
of the proper YM magnetic fields on a comoving scale L = 27r/fc in the position space is given by 



\B 
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k\C{k)\^ fc" 



1 
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27r2 



(2.12) 



where u^(= 1) and u'^{= 1) are the quantities denoting the dependence on the indices b and c, respectively. From 
Eq. ([TT2ll . we see that the YM magnetic fields evolves as \Bf''P'°P°'\t)\'^ = \B°-\'^/a'^, where \B°-\ is a constant. 

In this case, using the {fi, v) = (0, 0) component and the trace part of the (/x, v) = (i, j) component of Eq. 
where i and j run from 1 to 3, and eliminating I{R) from these equations, we obtain 
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(2.13) 



where H = d/a is the Hubble parameter and a dot denotes a time derivative, ' = d/dt. From Eq. (|2.8p . we see that 
the value of s depends on the field strength, in other words, it varies in time. The change in time of £, however, is 
smaller than that of other quantities because the dependence of e on the field strength is logarithmic, so that we can 
approximately regard s as constant in Eq. (j2.13p . In what follows, we regard e as constant. 
We consider the case in which f{R) is given by the Hu-Sawicki form [y, [^, 



fiR) = /Hs(i?) = 



ci {R/m^y 



C2 {R/iv?f + 1 



(2.14) 



which satisfies the conditions: limfl_nx) /hs(^) = Ci/c2 = const and limfl_,o fasiR) — 0. Here, ci and C2 arc dimen- 
sionless constants, n is a positive constant, and m denotes a mass scale. The above second condition means that 
there could exist a flat spacetime solution. Hence, because in the late time universe the value of the scalar curvature 
becomes zero, the YM coupling / becomes unity, so that the standard YM theory can be naturally recovered. To show 
that power-law inflation can be realized, we consider the case in which the scale factor is given by a{t) = d{t/tf, 
where t is some fiducial time during inflation, d is the value of a[t) at t — t, and p is a positive constant. In this 

case, H = p/t and R ~ 6 (h + 2H^^ = 6p(2p — l)/t^. At the inflationary stage, because R/m? ^ 1, we can use the 

. Substituting the above relations in terms of a, H and 

R, and the approximate expressions of /^g(i?), f^g{R) and f^g{R) derived from the above approximate expression 
of /hs(-R) for R/m? ^ 1 into Eq. (|2.13p . we find p = (?i + 1) /2. If n 1, p becomes much larger than unity, so 



approximate relation f^isiR) ~ (ci/c2) 1 — (1/^2) (R/m 



that power-law inflation can be realized. 
Eq. (|2.3p can be a source of inflation. 



Consequently, the YM field with a non-minimal gravitational coupling in 
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The constraint on a non-minimal gravitational coupling of matter from the observational data of the central tem- 
perature of the Sun has been proposed p^ . Furthermore, the existence of the non-minimal gravitational coupling 
of the electromagnetic field changes the value of the fine structure constant, i.e., the strength of the electromagnetic 
coupling. Hence, the deviation of the non-minimal electromagnetism from the ordinary Maxwell theory can be con- 
strained from the observations of radio and optical quasar absorption lines [sj , those of the anisotropy of the cosmic 
microwave background (CMB) radiation ISSl. 39l |. those of the absorption of the CMB radiation at 21 cm hypcrfine 
transition of the neutral atomic hvdrogen [4C | , and big bang nucleosynthesis (BBN) [U \^ as well as solar-system 
experiments [1^ (for a recent review, see j44l|). On the other hand, because the energy scale of the YM theory is higher 
than the electroweak scale, the existence of the non-minimal gravitational coupling of the YM field might influence 
on models of the grand unified theories (GUT). 



III. NON-MINIMAL VECTOR MODEL 



In this section, we investigate cosm olog y in non-abelian non-minimal vector model in the framework of F{R) gravity. 
The model action is given as follows [2l[ : 



Smg = d {Cmg + ^v) , (3.1) 



1 

1 



Cmg = 7^ {R + F{R)) , (3.2) 



Cy = I{R) i^--F;,F^^'' - V [A^']j , (3.3) 

where F{R) is an arbitrary function of R, F^^^ is given by Eq. (|2.4[) . and = g^^'^A'^A'^. {F{R) is the modified part 
of gravity, and hence F{R) is completely different from the non-minimal gravitational coupling of the YM field f{R) 
in jMl).) 

We should note that the last term V [A"^] in the action (|3.3|1 is not gauge invariant but can be rewritten in a gauge 
invariant way. For example, if the gauge group is a unitary group, we may introduce a (T-modcl like field U, which 
satisfies U^U — 1. The last term could be rewritten in the gauge invariant form: 

V [A"'^] V [-ctr {U^'D^U) {U^D^U)] , (3.4) 

where c is a constant for the normalization and Z?^ is a covariant derivative Z),j ~ d/^i + iA'lj^T'^ (T°'s are the generators 
of the gauge algebra). If we choose the unitary gauge U = 1, the term in (|3.4p reduces to the original one: V . 
This may tells that the action (|3.3p described the theory where the gauge group is spontaneously broken. 

Taking the variations of the action (|3.ip with respect to the metric f/^i^ and the vector field we obtain the 
gravitational field equation and the equation of motion of A'^ as 

(1 + F'{R)) R,,, - ^g^, {R + F{R)) + g^uUF'{R) - V^V,F'{R) = z^^yjv) ^ (3 5) 

CdV lA""^] 1 -\ 

+ i {f'{R)TR^,, + .g^.D [f\R)T] - V^V, [f\R)T] } , (3.6) 
.f = F^^i^''^'' + AV [A""^] , (3.7) 
I , , , ( , h dV\A'''^'^ ^ 



and ^=5m {^/^I{R)F''^"') - I{R) f^bc^b^pc^u _^ 2 ^ A"" = , (3.8) 



respectively, where T^if^ is the contribution to the energy- momentum tensor from A'^. 

We consider the case in which V is given by a class of the following power-law potential: V [A"^] = 

V {A^^/fh^)", where V is a constant, ni denotes a mass scale, and n{> 1) is a positive integer. Similarly to the 
preceding section, we consider the linearized equation of Eq. (|3.8p in terms of A° . For the above power-law potential, 
the form of the linearized equation of motion under the ansatz A'j{t, x) ~ and AQ{t, a;) = is the same as in the 
preceding section. (This is similar to the Coulomb gauge but since the action (|3.3p is not gauge invariant, or gauge 
symmetry is completely fixed by the unitary gauge as in after (j3.4p . this condition is only a working hypothesis.) 
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From = 0, we have (l/a^) AfAfdV [A^^] / = nV [A""^]. Using the = (0,0) component and the 

trace part of the (/x, ly) = component of Eq. (j3.5p . and eUminating I{R) from these equations, we obtain 



H + H^ + \ -F(R) - F'(R)H'^ + 3F"(R) 
' 6 



H + AlH^ + HH 



18F"'{R) (h + AHh\ 



f{R) (^-2H + 3i?2^ + 3/"(i?) [h - 2HH + AH^ - 2AH^H^ + 18 f"{R) {h + AHHy 



IB" 



+ 2 -/'(i?) 



nH + (1 + 2n - 2n^) H'^ 



3f"{R) H + 2(3-2n) HH + AH^ + 8 (1 - 2n) H^H 



18/"' (i?) (h + AHH) \v 



(3.9) 



In the case in which |Sf ^p^p'''^ (^p = \B''\'^/a^, V [A^^] cx a^^". If fi 2, the time evohition of V [A"'^] is the same 
as that of |Bf (P''°P°'''(<)|2. On the other hand, if n > 2, F [A"'^] decreases much more rapidly than |Sf (P™P'='')(i)|2 
during inflation. Hence, in the latter case we can neglect the terms proportional to V [A"^] on the right-hand side of 
Eq. (ifol) . 

We examine the following case. F{R) is given by [l3| 



F{R) 



[{R/AP) - (i?o/A/2)] 



2/+1 



{R^/M^ 



C3 + C4 { [(i?/M2) - (i?o/M2)]^'+i + {R^/Ap f+' } ' 



(3.10) 



which satisfies the conditions: limfl^oo F{R) = —M'^/ca = const, limij^o F{R) = 0. Here, C3 and C4 are dimensionless 
constants, I is a positive integer, and M denotes a mass scale. We consider that in the limit R — > 00, i.e., at the 
very early stage of the universe, F{R) becomes an effective cosmological constant, lim_R^oo F{R) = — AP/ca = — 2Ai, 
where Ai (^ ^^o^) is an effective cosmological constant in the very early universe, and that at the present time F{R) 



becomes a small constant, F{Ro) = — (^R^/M' 



, 21+1 



I c^+ca{Ro/M^) 



21+1 



= — 2i?o, where Rq {~ Ho ) is current 



curvature. Here, Hp is the Hubble constant at the present time: Hq = lOO/i km s ^ Mpc ^ = 2.1h x 10 ^^GeV 
1.5 X lO^^^eV [111, where we have used h = 0.70 [Ml. Moreover, /(i?) is given by [l3|: 



f{R) - ho{R) 



[{R/AP) (Ro/AP)] 



2q+l 



(Ro/M 



2N 2g+l 



C5 + C6 {[(i?/M2) - (i?o/M2)]2«+i + (i?o/Af 2)2^+1} 



(3.11) 



which satisfies the following conditions: lim^^tx) /no(-R) = 1/^6 = const, lim^^o /no(-R) = 0. Here, C5 and cg are 
dimensionless constants, and g is a positive integer. The form of F{R) in Eq. (|3.10[) and that of /no(-R) in Eq. (|3.11[) 
correspond to the extension of that of /hs {R) in Eq. (|2.14p . It has been shown in Rcf . [l3| that modified gravitational 
theories described by the action p.2|) with F{R) in Eq. p.lOp successfully pass the solar-system tests as well as 
cosmological bounds, and that they are free of instabilities. 

Using Eqs. ((3^ . ((3?T0)) and ((3TI|) . we find that at the very early stage of the universe {R/AP > 1 and R/AP > 

Rq/AP), a{t) cx exp ^^Ai/3t^ , so that exponential inflation can be realized, and that at the present time {F{R) = 
F{Ro) ~ — 2i?o), a{t) oc exp I y/Ro/3t] , so that the late-time acceleration of the universe can also be realized. 



IV. RECONSTRUCTION OF THE YM THEORY 



In this section, we indicate how to reconstruct the YM theory from the known evolution of the universe (for a 
review, see 



47[). We consider the following action: 



5 = 1 d'^V—gl^^+J^iF^^.F^'^-)^ . 

By introducing an auxiliary scalar field 0, we may rewrite the action ()4.ip in the following form: 



(4.1) 



(4.2) 
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Taking the variations of the action (|4.2p with respect to (p, we obtain 





dm pa pa,. I dQ{<t>) 



fJ.V 



(4.3) 



which could be solved with respect (f> as (j) = (j) {F^hF"" Here, the prime denotes differentiation with respect to 
Substituting the expression into the action (|4.2p . we obtain the action (|4.ip with 



p (i.a^F'^A.-) = - (P (0 (P;,F"'^'')) F^^.F^^'' + Q (0 . (4.4) 

Taking the variations of the action (j4.2p with respect to the metric tensor gi^i,, we obtain the Einstein equation: 

^ (r^, - l-Rg^J) = -\P{4>)F^^,FY + \g,u {P{^)F^,F^P^ + Q{4>)) ■ (4.5) 



On the other hand, taking the variations of the action (|4.2p with respect to , we obtain 

= a, (y^P((^)P'^''^) - ^Piq})r''^AlF^''^ . 



(4.6) 



For simplicity, we only consider the case in which the gauge algebra is SU{2), where /"■'"^ = e"^"^, and we assume that 
the gauge fields are given in the following form: 



.a _ J ac^(*)(5« (Ai = i = 1,2,3) 
^ " 1 (Ai = 0) 



(4.7) 



where a is a constant with mass dimension and A is a proper function of t. In general, if the vector field is condensed, 
the rotational invariance of the universe could be broken. In case of (|4.7p . the direction of the vector field is gauge 
variant. Hence, all the gauge invariant quantities given by (j4.7p do not break the rotational invariance. 
By the assumption, Eq. ()4.3p has the following form: 



a e a 



and {t, t)-component of Eq. (|4.5p is given by 



n lr2 3 ( -2{2„2\ -2 I -4 4A -4 

U = —^ri — -(aAe a + a c a 



The /X = component of Eq. (|4.6p becomes identity and fi = i component gives 



o = a 



(aP((/i)Ae 



3A 



(4.8) 



(4.9) 



(4.10) 



Since we can always take the scalar field (j) properly, we may identify the scalar field with the time coordinate (f) ~ t. 
By differentiating Eq. (|4.9p with respect to t and eliminating Q = dQ{<j>)/ (dt/)), we obtain 



= ^HH + a^X^e^^a^^P - P 



AA + A-' - X^H e^^a 



,2A„-2 



2a'' X-H] e^a 



,4A„-4 



(4.11) 



Furthermore, eliminating P by using Eq. (|4.10p . we find 



P = 



2HH 



2a2a-2e2^ A2 + AA 



Using Eq. (|4.12p . we can eliminate P (and P) in Eq. (|4.10p and obtain 

= 2 f A A + A^ + 'iX^'x] - a^e^^a-^H + 4 f A^ + AA - a^e^^Q-^H ] {X-H 



2 (A^ + AA) -a^e^^a-'^H 



2„2A„-2 1 



H H X ■ 

H H A 



2a^a 



2„-2„2A 



(4.12) 



(4.13) 
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If we give a proper a = a{t) and therefore H = H{t), Eq. (|4.13p can be regarded as a third order differential equation 
with respect to A. If we find the solution of A with three constants of the integration, we find the explicit form of 
P(0) = P{t) by using Eq. (|4.12p and further obtain Q(0) by using Eq. (|4.9p . Thus, we find the explicit form of three 
parameter families of the action (|4.2p . This tells that almost arbitrary time development of the university could be 
realized by the action (|4.2p or (|4.ip . 

As an example, we may consider the case of the power law expansion: a = [t/ti)^^ , H ~ hi/t, where ti and hi 
are constants. By assuming A = {hi — l)ln{t/ti) + Ai, where Ai is a constant, Eq. (14. 3p is reduced to the algebraic 
equation: = [2hi/ {hi - 1)] + {-Ahf + 13/ii + 2)X+{hi- if {hi - 2) (4i/i - 20), where X = aHfe^^. If this 
equation has a real positive solution with respect to X, we obtain Ai and therefore the exact form of A. Consequently, 
we can reconstruct a model to give the above power expansion. Similarly, any other evolutional history of the universe 
may be reproduced by the specific form of the action under consideration. 



FINITE-TIME FUTURE SINGULARITIES IN NON-MINIMAL MAXWELL-EINSTEIN THEORY 



In this section, we consider non-minimal Maxwell-Einstein theory with a general gravitational coupling. The model 
action is as follows 1301 : 



S'gr 



d'^Xy/^{CEli + Cem) , 



L 



1 



EM = -^I{R)F^ii^F'^'' , 



(5.1) 
(5.2) 



where Ceh is the Einstein-Hilbert action in 



F., 



d^A^ is the electromagnetic field-strength tensor. 



is the U{1) gauge field, and I{R) is an arbitrary function of R. 
Taking the variations of the action Eq. (|5.ip with respect to the metric g^i, and the U{1) gauge field A^, we obtain 
the gravitational field equation and the equation of motion of A^ as [22j 



R 



rp{EM) 



I{R) ( g'^^F^pF,^ - -g^.^F^pF'^" 



and 
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+ - I'{R)F^pF'^f'R^, + g^,n {l'{R)Fa,pF"P) - V^V„ {I' {R)F^pF"f) 
'^I{R)F^"') = , 



(5.3) 

(5.4) 
(5.5) 



respectively, where T^!^'^^ is the contribution to the energy-momentum tensor from the electromagnetic field. 

It follows from Eq. (|5.5p that in the flat FRW background, the equation of motion for the U{1) gauge field in the 
Coulomb gauge, Aj{t, x) = 0, and the case of Ao{t, x) = 0, is the same form as in Sec. II. We therefore obtain the 
approximate solution as \Ai{k,7])\^ = |C(fc)p, where |C(fc)p is given by Eq. (|2.1ip . Using this solution, we find that 
the amplitude of the proper magnetic fields in the position space is given by (t)p = (fc|C(fc)p/7r^) (fc''/a^) 

on a comoving scale L = 27r/fc. From this equation, we see that the proper magnetic fields evolve as (t)p ~ 

|i?P/a^, where \B\ is a constant. The conductivity of the universe Cc is negligibly small during infiation, because 
there are few charged particles at that time. After the reheating stage, a number of charged particles are produced, 
so that the conductivity immediately jumps to a large value: Uc 3> £[■ Consequently, for a large enough conductivity 
at the reheating stage, the proper magnetic fields behave in proportion to a~^(<) in the radiation-dominated stage 
and the subsequent matter-dominated stage . 

In this case, it follows from Eq. (|5.4p that the quantities corresponding to the effective energy density of the universe 
Poff and the effective pressure p^s are given by 



Pcff 



m 



+ H] I'{R) + QH [AHH + H] I"{R) 



6 



- 5H] I'{K) + 6 hOH^H - AH^ + HH - H] l"{R) 



36 [iHH 



(i?) 



(5.6) 
(5.7) 



We now suppose that I{R) is (almost) constant at the present time, and that for the small curvature, I{R) behaves 
as I{R) ~ IqR", where Iq and a are constants. We consider the case a < 0. The energy density of the magnetic fields 
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is given by ps = (1/2) |S/p™p'"') J(i?) = [\B\^ / (2a'')] /(i?). We take de Sitter background as the future universe. 
In such a case, when R becomes much smaller in the future, the energy density of the magnetic field becomes larger 
and larger in comparison with its current value. Thus, the strength of current magnetic fields of the universe may 
evolve to very large values in the future universe. 

In the flat FRW background, the Einstein equations are given by 3H^/k^ ~ pcs and — (^2H + 3H^^ j ~ p^ff, 

where peS and PcS arc given by Eqs. (j5.6p and (|5.7p . respectively. 

We examine the form oiI{R) producing the Big Rip singularity H ^ ho / {to — t) , where ho is a positive constant, and 
H diverges at t = to- In this case, the scalar curvature and the scale factor are given by i? ~ 6ho (2/io + 1) / (^o ~ i)^ 



and a ~ ao {to — t) ^° , respectively, where oq is a constant. We now assume that for the large curvature, I{R) 
behaves as I{R) ^ loR"'- Hence pcS in Eq. (j5.6p behaves as {to — t) 2"+""*"^ the left-hand side (l.h.s.) on the 
Fricdmann equation / = Pcs evolves as {to — t)^ ■ The consistency gives —2 = —2a + Aho, i.e., ho ~ {a — 1) /2 
or a = 1 + 2ho. The Friedmann equation 3H^/k'^ = pcs also shows = — [loho (l2ft.o + 6/io)" / (2ao), 

where we have used a = 1 + 2ho. This equation requires that /q should be negative. As a result, it follows from 
I{R) ~ /q-R" and a = l + 2ho that the Big Rip singularity can appear only when for the large curvature, I{R) behaves 
as If the form of I{R) is given by other terms, the Big Rip singularity cannot emerge. We note that if exactly 

J(i?) = IqR", H = ho/ {to — t) is an exact solution. 

Next, we investigate the form of I{R) giving a more general singularity H ~ ho{to—t)~^. In this case, 
the scalar curvature and the scale factor arc given by i? ^ Qho (3 + 2ho {to — t) {to — t) and a ^ 

'ho {p-i)-'{to-tr^^-'A 



respectively. 



ao exp 

In Ref. [49|; the finite-time future singularities has been classified in the following way: 

• Type I ("Big Rip") : For t ^ ts, a ^ oo, p ^ oo and \p\ oo. This also includes the case of p, p being finite 
at tg . 

• Type II ("sudden") [s^l : For t ^ tg, a ^ a^, p ps and \p\ oo 

• Type HI : For t ^ tg, a —^ a^, p ^ oo and \p\ — > oo 

• Type IV : For t ^ tg, a ^ Ug, p ^ 0, \p\ ^ and higher derivatives of H diverge. This also includes the case 
in which p {p) or both of p and p tend to some finite values, while higher derivatives of H diverge. 

Here, is, ag{y^ 0) and pg are constants. We now identify tg with to- The Type I corresponds to /3 > 1 or /3 = 1 case. 
Type II to — 1 < < case. Type HI to < /? < 1 case, and Type IV to /3 < — 1 but f3 is not any integer number. 
We note that if only higher derivatives of the Hubble rate diverge, then some combination of curvature invariants also 
diverges and it leads to singularity. 

We assume that for the large curvature, I{R) behaves as I{R) ~ loR"- If /3 < —1, in the limit t ^ to, R ^ 0- We 
therefore consider this case later, li (3 > 1, a oo and p^s and p^s — > because p^g oc a~^ and Peff ex. a~^. On 
the other hand, H oo. Hence the Einstein equations cannot be satisfied. 

If a > and < /3 < 1, peS in Eq. (|5.6p evolves as {to — i)^"''^^^\ but the l.h.s. on 3H^/k^ = pcs behaves as 
{to - ty'^'^. The consistency gives -2/3 = -a{l3+ 1), i.e., /3 = a/ (2 - a) or a = 2/3/ (/3 + 1). From "iH"^ / ^ pcff, 
we also find S/i^/k^ = - [/^ (6/io/3)" (1 - P) \B\'^] I [2a^ (/3 + 1)] , where we have used a = 2/3/ (/3 -f- 1) and on the l.h.s. 
we have taken only the leading term. This equation requires that /q should be negative. If a > and < /3 < 1, in 
the limit t ^ to, a ^ ao, R — *■ oo, poff — > oo and |pcff| — *■ oo. Hence the Type HI singularity appears. If a > and 
— 1 < /3 < 0, poff oo, but H ^ 0. Thus 3H^/k^ = pcs cannot be satisfied. 

If (/3 — 1) / (/3 -I- 1) < a < and — 1 < /3 < 0, in the hmit t ^ to, a ^ ao, R ^ oo, peS — > and |peff | oo. 
Although the final value of pcff is not finite but vanishes, this singularity can be considered to the Type II. The reason 
is as follows. In this case, when / and H are given by / = 1 + /oi?" and H = Ho + ho {to — t)~^ , respectively, where 
Ho is a constant, in the above limit pcff ~* Pa- From Eq. (|5.6p and 3H^ / = pcff, we find po — 3Hq/k^ = |Bp/ (2aQ). 
Hence po is a finite value. 

If a < {f3 — I) / {f3 + I) and — 1 < /3 < 0, in the limit t ^ to, a ao, R ^ oo, pcs ~^ and |poff | 0, but 

H ^ oo. Hence — cannot be satisfied. If a < and < /3 < 1, poff 0, but H oo. Thus 

3H^ / = PcS cannot be satisfied. 

In addition, we investigate the case in which f3 < —1. In this case, in the limit t to, a ^ ao and R ^ 0. We assume 
that for the small curvature, I{R) behaves as I{R) ~ /o-R"- If a > {/3 — I) / {P + 1), in the limit t to, Pcs 0, 
\pes\ — * 0, and higher derivatives of H diverge. Hence the Type IV singularity appears. If < a < {P — 1) / {j3 + 1), 

Pes and \pes\ oo. However, 7J — > and if — > 0. Thus — (2H + BiJ^j 1^ ~ PcS cannot be satisfied. 
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We remark that if I{R) is a constant (the case in which I{R) = 1 corresponds to the ordinary Maxwell theory), 
any singularity cannot appear. We also mention the case in which I{R) is given by the Hu-Sawicki form in Eq. (|2.14p 
as I{R) = fusiR) or expressed by Eq. ([3TT1l as I{R) = fm{R)- If /3 < -1 and I{R) = fusiR) or I{R) = ho{R), in 
the limit t to, a ao, i? — > 0, pcs — > and |pcff | ^ 0. In addition, higher derivatives of H diverge. Thus the Type 
IV singularity appears. 

As a consequence, it has been demonstrated that the Maxwell theory coupled non-minimally with the Einstein 
gravity may produce finite-time singularities in future, depending on the form of the non-minimal gravitational 
coupling. The general conditions for I{R) in order that the finite-time future singularities cannot appear are that in 
the hmit t ^ to, liR) ^ I, where /(^ 0) is a finite constant, I'{R) ^ 0, I"{R) and /'"(i?) ^ 0. 



VI. CONCLUSION 



In the present article, we have reviewed cosmology in non- minimal YM/Maxwell theory, in which the YM (elec- 
tromagnetic) field couples to a function of the scalar curvature. It has been shown that power-law inflation can be 
realized due to the non-minimal gravitational coupling of the YM field which may be caused by quantum correc- 
tions. Furthermore, we have considered non-minimal vector model in the framework of modified gravity. It has been 
demonstrated that both inflation and the late-time accelerated expansion of the universe can be realized. We have 
also discussed the cosmological reconstruction of the YM theory. In addition, we have studied late-time cosmology 
in the non-minimal Maxwell-Einstein theory. We have investigated the forms of the non-minimal gravitational cou- 
pling which generates the finite-time future singularities and the general conditions for this coupling in order that 
the finite-time future singularities cannot appear. Finally, it is interesting to mention that using the reconstruction 
method developed in Refs. [4^. [sij. one can present models with the crossing of the phantom divide in non- minimal 
YM-modified gravity generalizing the models suggested in Ref. [H] . 
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